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Abstract. Eight luminous quasars with —30 < Mg < 
—27 at z = 1.4 - 3.7 have been observed in the mid- 
and far-infrared using ISO. All the quasars have been de- 
tected in the mid-infrared bands of ISOCAM, while no 
far-infrared detections have been made with ISOPHOT. 
SEDs (Spectral Energy Distributions) from the UV to 
far-infrared have been obtained while supplementing ISO 
observations with photometry in the optical and near- 
infrared made from the ground within 17 months. The 
SEDs are compared with the MED (Mean spectral En- 
ergy Distributions) of low-redshift quasars with —27 < 
Mg < —22. It is shown that our far- infrared observations 
were limited by confusion noise due to crowded sources. 

Key words: Infrared:galaxies - Infrared:ISM:continuum 
- quasars: general - quasars: individual: PC 1548+437, 
PC 1640+4628, H 0055-259, UM 669, B 1422+231, PG 
1630+377, PG 1715+535, UM 678 - galaxies: photometry 



1. Introduction 

The IRAS sky survey was the first observation in which 
many quasars in the far-infrared were detected, and 
Neugebauer et al.(1986) presented IRAS measurements 
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of 179 quasars from 12 to 100 am. Combining IRAS ob- 
servations with data taken in other wavelengths, Sanders 
et al. ( |1989| ) have given SEDs from ~ 0.3 nm to 6 cm 
of 109 quasars in the Palomar-Green (PG) survey (Green 
et al. |1986[ ). Compiling Einstein, IUE, and IRAS data 
with a supplement of ground-based observations, Elvis et 
al. (1994) have presented SEDs for a sample of 47 nor- 



mal quasars, and derived the SEDs for radio-quiet and 
radio-loud quasars. Recently, Andreani et al.( 1999Q have 
presented SEDs for 120 optically selected quasars at low- 
and high-redshift, combining sub-mm and mm observa- 
tions with optical, near-infrared, IRAS, and radio obser- 
vations. 

The gross shape of quasars SEDs is characterized by 
two major features; the big blue bump shortward of 0.3 fim 
and the infrared bump between 2 and 200 /im (Sanders et 
al. 1989 ; Elvis et al. 1994 ). It is generally considered that 
the blue bump is dominated by thermal emission from an 
accretion disk. The infrared bump is ubiquitous in both 
radio-quiet and radio-loud quasars. The infrared emission 
in radio-loud quasars has generally been taken to be domi- 
nated by non-thermal emission (e.g., Impey & Neugebauer 



198$ Bloom et al. |1994| ; Neugebauer & Matthews |1999J). 



The origin of the infrared emission in radio quiet 
quasars is generally attributed to thermal emission from 
heated dust. The rise from 1 um minimum toward 3 um, 
which is universally present in quasars, is naturally ex- 
plained by the sublimation of dust grains at ~ 1,500 K 
(e.g., Kobayashi et al. 1993). This is expected from the 



current unified models of active galactic nuclei, which have 
a dusty obscuring torus around a central source. Pier & 



Krolik (1993) computed the infrared properties predicted 



for the dust obscuring tori surrounding central sources 
by modeling free parameters of the inner radius of the 
torus to its thickness, the Thomson depths constraining 
the outer radius of the torus, and the flux of the nu- 
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clear radiation. The infrared emission in 1 - 10 /jm of 
PG quasars is well explained by their models, while the 
emission in the far-infrared is much greater than that pre- 
dicted. Several models were proposed to account for the 
far- infrared emission of quasars; these are (l)warm dust in 
a distorted disk extending from O.lkpc to more than lkpc, 
which is heated directly by radiation from central source 
(Sanders et al. 1989| ), (2)dust clouds in the narrow line re- 
gion heated by the radiation from the centra l sources and 
star- forming regions (Rowan- Robinson 1995 ), and (3)dust 
in the obsc uring tori extending to 200 - 3kpc (Andreani 
et al. 1999] ). However, it should be noted that temporal 



variations of 10 fxm brightness of radio-quiet quasar PG 
1535+547 found by Neugebauer & Matthews ( |1999[ ) casts 
doubt on the thermal origin of the infrared emission in 
radio-quiet quasars. 

High-redshift quasars have been widely recognized to 
provide unique probes of high-redshift star formation 
and galaxy evolution. Metallic abundance in high-redshift 
quasars is solar or higher metallicities out to z > 4. Com- 
paring the UV-to-optical spectra of 186 quasars with < 
z < 3.8, Osmer et al. ( [1994 ) concluded that there was no 
evidence for redshift-dependent spectral changes. Hamann 
k Ferland Ql993j ) analyzed N V /C iv and N V /He n broad 
emission ratios, and found that metallicities in the broad 
line gas at high redshift are 1-10 times solar. The ra- 
tio of UV Fe n/Mg ii of B 1422+231 at z = 3.6 shows the 
the host galaxy was already in the late-evolutionary phase 
of th e Fe enrichm ent by SNe la at z = 3.6 (Kawara et al. 



1996; Yoshii et al. 1998). These all suggest that the central 



part of host galaxies formed rapidly at very high-redshift 
z > 10, and metallicity enrichment in the central part was 
already completed at z = 4 — 5. 

Dust emission from high-redshift quasars provide an- 
other probe to study star formation and evolution of host 
galaxies, especially in outer regions. Up to date, many 
quasars at z > 3 have been detected in the region from 
350 fim to 1.3 mm (e.g., Andreani et al. 1993|; Chini & 
Kriigel |1994[ Du nlop et al. |1994| ; Isaak et a l. |1994| ; McMa- 
hon e t al. |1994[ Ivison |1995|; Om ont et al. [19964 Hu ghes 



et al. |1997[ Benford et al. |199J; Carilli et al. p000|) . An 



extensive 240 GHz survey by Chini et al. ( 1989|) revealed 
that the majority of z < 1 quasars have dust masses about 
a few times 10 7 M Q , comparable to normal spiral galaxies, 
thus suggesting that dust is heated by radiation from the 
central sources. On the other hand, dust masses > 10 8 M© 
have been found in six of 16 quasars at z > 4 (Omont et 
al. |!996a| ), and CO emission was found in three of these 



(Ohta et al. 1996; Omont et al. 1996b Guilloteau et al 



1997,1999; Carilli et al. 1999). This may imply that dust 



emission in these host galaxies may be dominated by ra- 
diation from star- forming regions at high-redshift. 

To study the distribution of dust in high-redshift 
quasars, mid- and far-infrared observations are indispens- 
able. The emission from sublimated dust could be ob- 
served in the mid-infrared, thus probing that the inner- 



most part of the obscuring tori, and the distribution of 
dust from obscuring tori to disks (or outer star-forming 
regions) could be obtained from far-infrared observations. 
The IRAS sky survey has only detected several quasars 
at z > 3 with exceptionally high luminosity in the far- 



infrared (Neugebauer et al. 1986; Bechtold et al. 1994; 
Irwin et al. |l~998| ); F08279+5255 (APM 08279+5255) at z 
= 3.87, 2126-158 (PKS 2126-15) at z = 3.275, and 0320- 
388 at z = 3.12. In hoping to detect more quasars, several 
groups have carried out mid- and far-infrared photome- 
try us ing t he Infrared Space Observatory (ISO; Kessler 
et al. 1996 ). The first report of the ISO European Cen- 
tral Quasar Programme which observed 70 quasars be- 
tween 4.8 and 200 [im includi ng hig h-redshift quasars has 
been published by Haas et al. ( |1998|) , and the I SO /NASA 
AGN Key Project has also observed 72 quasars and AGN s 
covering a range of redshift up to 4.7 (Wilkes et al. |1999| ). 
As an ISO open time program, we have executed mid- and 
far-infrared observations of eight quasars at 1.4 < z < 3.7 
using the raster mapping mode. This paper presents the 
results of the photometry of these quasars supplemented 
with optical and near-infrared data taken on the ground. 
In the course of this work, mid- and far-infrared sources 
have scrcndipitously been discovered. These sources will 
be described in the forthcoming paper (Oyabu & Kawara 
|2000p . Throughout this paper, H = 75km s" 1 Mpc- 1 
with qg = is assumed. 

2. Observations and Results 

Eight quasars have been observed with ISO in revolu- 
tion 169 - 781 (1996 March - 1998 January). Optical and 
near-infrared imaging observations have followed within 
24 months (mostly within 17 months). 

2.1. The Sample 

Table [l] lists the sample of the eight quasars in the 
sequence of ISO revolutions for execution. The sample 
consists of luminous quasars with Mb < -28 except PC 
1548+4637 and PC 1640+4628; these two quasars which 
were observed at the beginning of this work, turned out 
to be too faint to be detected in the far-infrared, and thus 
the sample selection criterion was changed to include 
very luminous quasars in low far-infrared background of 
infrared cirrus emission. Fig. |l| presents the eight quasars 
(large filled diamonds) on the z — Mb plane together 
with t hose (small filled circles) in the sample by Elvis 
et al.( 1994[ ) and those ( faint gray points) complied by 
Veron-Cetty & Veron ( 1998 ). All the sample quasars 
are radio-quiet or optically selected except B 1422+231 
which is a core-dominant flat-spectrum radio source. 

2.2. Mid-infrared observations with ISOCAM 

The mid-infrared observati ons w ere performed with 
ISOCAM (Cesarsky et al. |l996|) . Three broad band 
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Table 1. Quasars' Sample 



Object 


R.A. (J2000) Dec. 


z 


M B a 


Radio" 


Rev (UT) C 


Other name and Notes 


PC 1548+4637 


15:50:07.6 


+46:28:55 


3.544 


-27.0 


Quiet 


169 (960504) 




PC 1640+4628 


16:42:05.1 


+46:22:27 


3.700 


-26.8 


Quiet 


185 (960520) 




H 0055-2659 


00:57:58.1 


-26:43:14 


3.662 


-29.2 


Optical 


380 (961130) 




UM 669 


01:05:16.8 


-18:46:42 


3.037 


-28.4 


Optical 


415 (970104) 


Q 0102-190 


B 1422+231 


14:24:38.1 


+22:56:01 


3.62 


-29. 8 d 


Loud 


424 (970113) 


Lensed quasar 


PC 1630+377 


16:32:01.1 


+37:37:49 


1.478 


-28.2 


Quiet 


424 (970113) 


Also observed on rev. 778 (980101) 


PC 1715+535 


17:16:35.4 


+53:28:15 


1.940 


-28.5 


Quiet 


712 (971027) 




UM 678 


02:51:40.4 


-22:00:27 


3.205 


-29.4 


Optical 


781 (980104) 


Q 0249-222 



The absolute B magnitude for Ho = 75km s~ Mpc~ with qo = 0.0. 

Radio property; Quiet = radio quiet, Loud = radio loud, and Optical = optically selected. 
UT is given in the yymmdd format where yy = ye ar, m m = month, and dd = day. 
Mb ~ ~26 after demagnification (Kormann et al. 



1994) 



Table 2. Filters and settings for raster mapping 



Filter 2 


LW2 


LW3 


LW10 C_90 


C.160 


Are/ 


6.7/im 


14.3/xm 


12.0/im 90^m 


170 /j,m 


AA 


3.5/im 


6.0/im 


7.0/im 51/im 


89/rni 


Airy diameter 6 


5.6" 


12.0" 


10.1" 76" 


143" 




Raster points M x N, Raster step, Exposure time per raster point c , Pixel field of view 


QSO field 


LW2 


LW3 


LW10 C_90 


C.160 


PC 1548+4637 


5x2, 6", 60s, 6" 




4x1, 6", 25s, 6" 4x4, 44", 64s, 44" 


3x3, 90", 64s, 89" 


PC 1640+4628 


5x2, 6", 60s, 6" 




4x1, 6", 25s, 6" 4x4, 44", 64s, 44" 


3x3, 90", 64s, 89" 


H 0055-2659 


5x2, 6", 60s, 6" 




5x3, 44", 69s, 44" 


4x2, 90", 73s, 89" 


UM 669 


6x6, 7", 90s, 3" 






10 x 10, 46", 16s, 89" 


B 1422+231 


4x4, 7", 90s, 3" 


4x4, 7", 90s, 3" 




10 x 10, 46", 16s, 89" 


PC 1630+377 d 


4x4, 7", 38s, 3" 


4x4, 7", 38s, 3" 




10 x 10, 46", 16s, 89" 


PC 1715+535 


4x4, 7", 38s, 3" 


4x4, 7", 38s, 3" 




10 x 10, 46", 16s, 89" 


UM 678 


4x4, 7", 90s, 3" 


4x4, 7", 90s, 3" 




6x5, 92", 64s, 89" 



for 



Cited from ISO Handbook Volume III (CAM) (Siebenmorgen et al. |l99J) and Volume V (PHT). 
The aperture photometry was performed by using the apertures with a diameter of 2 x dairy for ISOCAM and of dairy 
ISOPHOT. The diameter of the Airy disk d alr y is 2.44(A/60cm) in radian or 0.84A(/xm) in arcsec. Note that the FWHM of 
the Airy disk is 0.422 x dairy 

ISOCAM exposure time per raster point is given as Ti nt x N exp where Ti n t is an integration time of a single exposure and 
N exp is the number of exposures. Ti„t = 5sec was used for all LW2 and LW3 observations except for PC 1630+377 and PC 
1715+535 for which Ti„ t = 2sec was used. Tint = 2sec was used for LW10 observations. 

The second observation was executed on revolution 778 in C_160 with a parameter set of (6 x 6, 92", 64s, 89"). 



filters, namely, LW2 (reference wavelength 6.7/mi), 
LW3(14.3yum), and LW10 (12.0/im) were used. All the 
quasars were observed in LW2 with additional measure- 
ments in LW3 or LW10. To detect faint sources against the 
intense background dominated by zodiacal light, the AOT 
(Astronomical Observation Template) CAM01 which is 
the microscan raster mapping mode was used to achieve 
accurate flat-fielding (e.g., Taniguchi et al. 1997; Altieri et 
al. |1999[ ). Table | lists the details of CAM01 parameters 
as well as the characteristics of the broad band filters. 

The standard ISOCAM reduction software CIA 3.0f] 
was used to produce ISOCAM images from ERD (Edited 



1 CIA is distributed by the ISO Data Centre and the ISO- 
CAM Data Centre at Saclay, France. 



Raw Data). This process includes dark subtraction, 
deglitching, correction for the transient behavior of ISO- 
CAM pixel signals, and flat fielding(Delaney |1998| ) The 
inversion transient correction model of Starck et al.(1999) 
was applied. The factor of the correction for the tran- 
sient behavior is 0.58 - 0.77 in LW2, 0.79 - 1.0 in LW10, 
and 0.87 in LW10. Fig. | shows LW2 and LW3 maps for 
the two brightest quasars B 1422+231 and PG 1715+535. 
All the quasars were clearly detected at the expected po- 
sitions. Aperture photometry was performed using IDL. 
Two apertures centered on the object were used; the small 
one has a diameter of 2d a iry, two times the Airy diame- 
ter as given in Table and the other has a diameter of 



4d„ 



The photometry was corrected for loss of flux in 



the PSF (Point Spread Function) wings by computing the 
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Table 3. ISO mid/far- infrared flux density and UV to IR luminosity 







ISOCAM(mJy) 




ISOPHOT(mJy) 


Luminosity(10 


13 r \a 


Object 


LW2 


LW3 


LW10 


C_90 


C.160 


Luvo 




L\.2hp.m 


PC 1548+4637 


0.16 + 0.04 




0.67 + 0.17 


8.2 ± 28 


32 + 80 


11 


< 130 


1.7 


PC 1640+4628 


0.15 + 0.04 




0.29 + 0.11 


-8.3 + 20 


83+121 


8.1 


< 240 


1.8 


H 0055-2659 


0.29 + 0.06 






16 + 19 


-2.7 + 60 


24 


< 160 


3.3 


UM 669 


0.50 + 0.01 








-9.9 + 18 


27 


< 39 


3.5 


B 1422+231 


5.8 + 0.06 


15.1 + 0.07 






83 + 56 


410 


< 320 


64 


PC 1630+377 


4.3 ± 0.04 


7.3 + 0.07 






5.6 + 14 6 


28 


< 16 


4.3 


PC 1715+535 


4.7 + 0.06 


10.3 + 0.06 






-13 + 16 


81 


< 26 


9.4 


UM 678 


0.73 + 0.01 


1.9 + 0.03 






-15 + 58 


30 


< 116 


5.9 



a The UVO (0.1/im to 1/im) luminosity Luvo = 47rD| x 2.3 J^^uf u d\og l0 u, the IR (1.5/im to lOOfim) luminosity Li r = 

o rl4 3 o 

A-kDI x 2.3 J 12 S WUlog 10 v, and the 1.25 (im luminosity L\.2z^m = iwD L uf l ,, where Dl and /„ denote the luminosity 
distance and observed flux density, respectively. Where no observations exit such as 1/im flux denisities, interpolated flux 
densities were used. Luvo is based on the optical and near-infrared flux densities given in Table ^. 
6 /i/(174^m) = 89 + 78 mjy was obtained in the second observation. 



-32 



2.3. Far-infrared observations with ISOPHOT 



The far-infrared observations were made with ISOPHOT 



-30 - 



-28 



-26 



-24 



-22 



Fig. 1. Sample quasars (large filled diamonds) compared 
with other samples of quasars on the z - Mb plane. Small 
filled circles show quasars studied by Elvis e t al.(1994 ). 
Quasars cataloged by Veron-Cetty & Veron ( 1998 ) were 
plotted by dots, which appears as a gray background on 
this plane. 



PSF based on the model having a two mirror f/15 tele- 
scope with radii for the primary and secondary mirror of 
30 and 10 cm, respectively(Muller |l999| ). The factor of 
the PSF correction is 0.75 - 0.90 (i.e., loss of flux is 0.1 - 
0.25), depending on the raster step and the pixel field of 
view. 



The results after these corrections are given in Table 
H with statistical errors. The errors in the absolute photo- 
metric calibration are not included in Table 3; these errors 



are estimated to be 15%(Siebenmorgen et al. 1999). 



(Lemke et al. 1996). All the quasars were observed in the 
broad C_160 (170/im) band with additional measurement 
in C_90 (90/im). ISO far-infrared surveys (Kawara et al. 
1998; Puget et al. 1999| ) clearly indicate that the sky seen 



in the far-infrared has a clumpy structure which is made 
up of IR cirrus and extragalactic sources. This structure 
rotates with time relative to the ISO coordinate system 
due to the field rotation, increasing the probability of fault 
detection if the chopping mode is used. We thus selected 
the PHT22 staring raster map mode to make small maps 
around the quasar. Table ^ presents the characteristics 
of the photometric filters and the details of observational 
parameters for raster mapping. 

ISOPHOT images were produced by using the stan- 
dard ISOPHOT redu ction software PIA V7.3 and V8.l| 



(Gabriel et al. 1997), starting at the edited raw data 



(ERD) created via the off-line processing version 7.0. The 
AOT/Batch processing mode of PIA is used with the 
default parameters to reduce ERD to the Astronomical 
Analysis Processing (AAP) level. This standard reduction 
includes linearization and deglitching of integration ramps 
on the ERD level, signal deglitching and drift recognition 
on the SRD (Signal per Ramp Data) level, reset inter- 
val normalization, signal deglitching, dark current sub- 
traction, signal linearization, and vignetting correction^ 
on the SCP (Signal per Chopper Plateau data) level. The 
responsivity calibration was made on the SPD (Standard 
Processed Data) level by using the second measurement of 



2 PIA is a joint development by the ESA Astrophysics Divi- 
sion and the ISOPHOT consortium. 

3 Vignetting correction is set to on in the default setting of 
PIA. However, the effect of the correction is none; just multi- 
plying the data by 1.0, because the chopper throw is zero in 
the PHT22 raster mode. 
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1422+231 



PG 1715+535 



O 
CO 




Fig. 2. ISOCAM LW2(6.7^m) and LW3(14.3/xm) and 
ISOPHOT C_160(170/mi) maps are shown on the top, 
middle, and bottom panels, respectively. B 1422+231 
shown on the left was observed on revolution 424 and 
PG 1715+535 on revolution on 712. The ISOCAM and 
ISOPHOT maps are 90"and 598"(10'), respectively, with 
North and East as indicated by the arrows. Quasars should 
be at the center of the gray circles. 



the internal Fine Calibration Source 1(FCS1) which is cal- 
ibrated against celestial standards. The correction for the 
transient behavior of the detectors was applied to point 
sources (quasars) on this level. Images were produced on 
the AAP (Astronomical Analysis Processing) in the map- 
ping mode with median brightness values. 

The correction for drift in the responsivity is not im- 
portant to our observations, and so this correction was not 
made. To check the importance of the drift, the MEDIAN 
filter technique was applied to the results from AAP (here- 
after called AAP map). Applying this technique to large 
AAP maps in the Lockman hole, Kawara et al. (1998) 
show that this is a powerful tool to correct for drift in the 
detector responsivity. However, unlike large AAP maps in 
the Lockman hole, the MEDIAN filter technique does not 
improve our AAP maps. This is attributed to the size of 
our AAP maps; the observing time of these small maps is 



shorter than the timescale of drift in the detector respon- 
sibility, and so impact by the drift is small. In fact, every 
detector pixel of the C100 and C200 detector arrays has 
been checked for signal, and neither spike noise nor drift 
in the responsibility was found. Maps with all the detec- 
tor pixels always give better results than those obtained 
by masking some of the detector pixels. In addition, it 
was confirmed that there were no significant differences 
between two AAP maps produced by two different algo- 
rithms on the AAP level, namely, the full coverage and 
distance weighting algorithms. Fig. || shows AAP maps 
of C_160 for the brightest quasars B 1422+231 and PG 
1715+535. 

Aperture photometry was then performed using 
IRAFg and Skyviewf] in the manner similar to ISOCAM. 
The original pixel sizes of AAP maps are equal to the 
raster steps. Because these are too large to center the aper- 
ture on the quasar with sufficient accuracy, AAP maps 
were rebinned in such a way that each original pixel is 
converted into 10 x 10 sub-pixels. Two apertures centered 
on the object were used; the small one has a diameter of 
dairy, the Airy diameter as given in Table ||, and the other 
has a diameter of 2d a i ry . The photometry was corrected 
for loss of flux in the PSF (Point Spread Function) wings 
by computing the same PSF model as used for ISOCAM. 
The factor of the PSF correction is 0.63 except for PC 
1548+4632 and PC 1640+4628. These two quasars were 
centered between four pixels in such a way that quasars 
illuminate the four pixels equally. Consequently the loss 
of flux measured with the two apertures is large, and the 
factor of the PSF correction is 0.27 for C_90 and 0.23 for 
C.160. 

The results after these corrections are given in Table 
H with statistical errors. The errors in the absolute photo- 
metric calibration are not included in Tabled; these errors 
are estimated to be 30% (Klaas et al. 2000] ). It is noted 
that PG 1630+377 was observed twice at 170^m to check 
variability. 



2.4- Optical and near-infrared data from ground-based 
observations 

Optical images were taken on the 0.9m telescopes at CTIO 
and the Schmidt 1.05m telescope at Kiso Observatory. 
Near-infrared imaging was made in the standard dithering 
mode on the 1.88m telescope at the Okayama Astrophys- 
ical Observatory, NAOJ. SExtractor (Bertin & Arnouts 
1996 ) was applied to the optical and near-infrared images 
to perform photometry. Flux cal ibrati on was made using 
the standards given by Landolt (|l992|) in the optical and 
the UKIRT faint standards (Casali et al. |1992 ) in the near- 
infrared. A typical photometric error is 0.05 mag. 



4 IRAF is distributed by NO AO. 

5 Sky view is distributed by IPAC. 
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Table |J presents magnitudes of the quasars together 
with statistical errors. As shown in the table, all the 
ground-based observations were performed within 24 
months (mostly 17 months) from the ISO observations so 
as to reduce the chance of having flux variations between 
ISO and ground-based observations. Table |] also supple- 
ments optical magnitudes taken by Richards et al. (1997) 
and Schneider et al. (1994). All the supplementary quasars 
but PG 1548+4637 were observed within 27 months before 
the ISO observations. 



3. Discussion 



All the quasars have been detected in the ISOCAM LW2 
and LW3 or LW10 band. Combining optical and near- 
infrared observations, SEDs from the UV to mid-infrared 
have been obtained for eight quasars at z = 1.4 - 3.7. Fig. 
U shows the restframe SEDs (£„[J. The down arrows in- 
dicate the 3 a upper limits obtained with the ISOPHOT 
observations. Note that 3 a upper limits in Fig. |^ cor- 
respond to three times uncertainties listed in Table |j| 
The mean spectral distributions (MEDs) of low redshift 
quasars by Elvis et al.(1994) are also shown by fitting to 
the ISO 6.7/im flux density. Note that 6.7/im in the ob- 
serving frame is approximately 1.25/im in the restframe. 
The dash lines indicate the 1 a deviation (strictly speak- 
ing 68 Kaplan-Meier percentile) envelopes. The rise from 
the 1 /im minimum toward 3 /zm is evident at least in 
z = 3.5 quasar PC 1548+4637. This rise is naturally ex- 
plained by the sublim ation of dust grains at ~ 1,500 K 
(e.g., Kobayashi et al. 1993 ) as expected from the current 
unified models of active galactic nuclei. Hence, our data 
suggest that the obscuring torus model could be applied 
to high-rcdshift quasars. 

The original objective of this work was to study 
whether there was any evidence for redshift-dependent 
changes in SEDs from the UV to far-infrared. As shown in 
Fig. the data from the UV to mid-infrared do not devi- 
ate significantly from the MED of the low-redshift sample. 
However, this cannot be taken as evidence for no SED evo- 
lution in the UV to mid-infrared, because our sample is 
too small in number and too narrow in spectral coverage. 
Our far-infrared observations only provide 3 a upper lim- 
its, and are not very helpful to check SED evolution in the 
far-infrared. 

The far-infrared detection with ISO have been 
severely limited by the confusion due t o gala xies and local 
peaks of th e IR cirrus (Ka wara e t al. 1998 ; Herbstmeier 
et al. [199$ Matsuhara et al. [20001 ). Simulating the 90 and 
170 /im source counts in the Lockman Hole, Kawara et al. 
( 2000a , 2000b ) concluded that the effect of confusion due 
to crowded sources (presumably extragalactic) was signifi- 
cant in the flux range below 200 mJy, and 50% of 140 mJy 



sources at 170 jum were left undetected due to confusion 
noise. The expected fluxes of our sample quasars ranges 
from 7 - 170 mJy at 170 /im. Such confusion is clearly 
seen in Fig. ||. There are a few sources scattering around 
the quasar, which hampers detection in the far-infrared. 
It should be noted that the brightest source in the PG 
1630+377 field has 106 mJy and the o ne in PG 1715+535 
has 190 mJy (Oyabu & Kawara |2000|) . 



It is clear that larger aperture telescopes like FIRST, 
which provide spatial resolution finer than ISO, are re- 
quired for far-infrared studies of high-rcdshift quasars. 
In cases where similar aperture telescopes such as 
IRIS(ASTRO-F) and SIRTF are used, data sampling, 
namely a step size of raster mapping, should be care- 
fully designed so that far-infrared images can be improved 
by deconvolution without violating the sampling theorem 



(Magain, Courbin, & Sohy 200C) 



6 L v is defined by A-kD\v] v where Dl denotes the luminosity 
distance to the object. 
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Universal time when observed in the yymmdd format. 

At the Kiso Observatory, the 2048 2 CCD with 1.5"per pixel was used on the 105cm Schmidt telescope. 
At CTIO, the 2048 2 CCD with 0.4" per pixel was used on the 90cm telescope. 

At OAO (Okayama Astrophysical Observatory), the infrared imager spectrometer which has a HgCdTe 256 2 detector array 
with 0.97"per pixel was used on the 188cm telescope. 

All but PC 1548+46 37 we re observed within 27 month s before the ISO observations. 
R = Richards et al. ( {L997[ ); S = Schneider et al. ( [19941 ). 
The T4 bandpass was used. 
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Fig. 3. Rcstframe SEDs compared with the MED of low-rcdshift quasars compiled by Elvis et al. ( 1994 ). Dashed lines 
show the 1 a deviation (actually 68 Kaplan-Meier percentile) envelopes. Asterisks denote data taken by Richards et 
al.( [L997D and Schneider et al .( |l994| ). The dot-dashed and dot lines f or B 1422+231 shows the MED of radio-loud 
quasars from Elvis et al.(1994) and the SED of 3C273 from Lichti et al.(1995), respectively. Note that 3 a upper limits 
correspond three times uncertainties listed in Table [| 
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